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ABSTRACT 

We examine the scattering of acoustic p-mode waves from a thin magnetic 
fibril embedded in a gravitationally stratified atmosphere. The scattering is me- 
diated through the excitation of slow sausage waves on the magnetic tube, and 
only the scattering of the monopole component of the wavefield is considered. 
Since such tube waves are not confined by the acoustic cavity and may freely 
propagate along the field lines removing energy from the acoustic wavefield, the 
excitation of fibril oscillations is a source of acoustic wave absorption as well as 
scattering. We compute the mode mixing that is achieved and the absorption 
coefficients and phase shifts. We find that for thin tubes the mode mixing is weak 
and the absorption coefficient is small and is a smooth function of frequency over 
the physically relevant band of observed frequencies. The prominent absorption 
resonances seen in previous studies of unstratified tubes are absent. Despite the 
relatively small absorption, the phase shift induced can be surprisingly large, 
reaching values as high as 15° for / modes. Further, the phase shift can be 
positive or negative depending on the incident mode order and the frequency. 

Subject headings: MHD — Sun: helioseismology — Sun: magnetic fields — Sun: 
oscillations 



1. Introduction 



It is well established that sunspots and magnetic plage are strong absorbers and scatter- 
ers of the sun's acoustic p-mode waves (e.g., lBraunlll995l ; iBraun fc Birchl 120081 : lHaber et al. 
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19991 ). The absorption is unexpectedly large, reaching values as high as 70% for sunspot 



groups and a lesser, but still substantial, value of 20% for magnetic plage. Equally strong 
phase shifts are also observed for sunspots, with shifts as large at 100° being possible. Inter- 
estin gly, plages lack measurable phase shifts, despite possessing significant absorption (IBraun 



1995|). 



Initial attempts to model the interaction of acoustic waves with magnetic features con- 
centrated on the absorption caused by sunspots, as the observed signal was strong and easy 



to measure. A host of a 
tion (e.g 



p sorption mechanisms were inves tigated, includin g resonant ab sorp- 



Hollwed 11988 



Loul 



1990 



Keppens et al 



mode conversion (e.g., ICally &: Bogdanlll993l ; ICally 



99I, mode mixing (ID'Silval Il994l ) and 
2000). The present belief is that mode 



conversion from acoustic wav es to slow magnetoson ic waves within an appreciably inclined 
magnetic field is responsible (jCrouch fc Callyl 120031 ) . Under this paradigm, acoustic energy 
is transformed into slow waves which propagate along the field lines removing energy from 
the p-mode cavity. The conversion process is most efficient when the phase velocity of the 
incident acoustic wave matches the local phase velocity of the slow mode. Such strong cou- 
pling of the phase velocities occurs when two conditions are met: the wavenumber vector of 
the incident acoustic wave is nearly parallel to the magnetic field and the Alfven and sound 
speeds are comparable. In a sunspot, the Alfven and sound speeds are equal to each other 
in a layer just below the solar surface and the ray paths of the p modes in this layer are such 
that the strongest conversion occurs in regions with field lines that are inclined 20-30 degrees 
from vertical (such as within penumbrae). In addition to successfully explaining the absorp- 
tion observed for s unspots, mode co nversion models have also succeeded in reproducing the 
large phase shifts ( ICally et al.l 120031 ) . 

So, while we presently have a fairly good understanding of how a flaring, monolithic, 
magnetic structure such as a sunspot can absorb and scatter acoustic waves, we have made 
little progress in modelling the observed absorption by plage. The absorption within a plage 
is smaller than that of a sunspot by only a factor of two, yet the mode conversion process 
is unlikely to work in the same manner, as the field within a plage has a fractured and fibril 
nature, and probably has less inclined field of a structured geometry in the photospheric lay- 
ers. The limited theoretical work that has been done modelling the effect s of fibril magnetic 
fields indicates that multiple scattering may pla y a very important role ( IBogdan &: Zweibel 



19871 : Keppens et al.lll994i IBogdan fc Foxlll99ll ). iKeppens et al.l ( 119941 ) examined multiple 



scattering between a small number of flux tubes (2-3) in an unstratified atmosphere and 
found that resonant abso rption by th e collection of tubes can be greatly enhanced if the 
tubes are tightly spaced. iTirryl (120001 ) performed a similar calculation but considered inci- 
dent p-mode waves trapped within a solar-like wave guide and their interaction with thin 
flux tubes. Tirry permitted tube waves excited on the magnetic fibrils to freely cross the 
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bottom boundary of the acoustic waveg uide and discovered th at such wave leakage destroyed 
the absorption resonances observed by iKeppens et al.l (119941 ) . Both studies lacked gravita- 
tional stratification. Neither of these studies included a sufficiently large number of tubes to 
reliably estimate the expected phase shifts measured in plage. The lack of measurable phase 
shift is likely to be the result of incoherent scattering from a large number of tubes. 



Jain et al.l ( 120091 ) examined whether the wave leakage mechanism explored by Tirry is 



sufficiently large on its own to explain the absorption observed within plage. They modeled 
a plage as a large collection of thin magnetic tubes embedded in a stratified atmosphere 
and considered each tube to absorb acoustic waves in isolation from the other members of 
the collective. The MHD tube waves excited on each tube by the incident p mode freely 
propagate along the fibril out of the acoustic cavity despite the gravitational stratification. 
They found that high levels of absorption (potentially greater than 50%) were easily achieved. 
However, they did not take into account scattering between tubes or the screening that is 
likely to occur for tubes located far from the edge of the plage. 

To date few studies have examined both the absorption and scattering of acoustic waves 
by magnetic tubes res iding in a gravitationally stra tified atmsophere. Notable are those of 
Hanasoge et al.l (120081 ) and lHanasoge k, Callyl (120091 ) which focused on the scattering of the 
dipole component of incident f modes from thin mag netic tubes through the excitation of 
kink waves, and the study by iGordovskyy et al.l (120091 ) that, through numerical simulation, 
examined axisymmetric scattering from tubes of general thickness. Here we examine scat- 
tering and absorption of the wavefield's monopole component of both / and p modes from 
a thin magnetic tube in a gravitationally stratified atmosphere. The scattering is mediated 
through the excitation of axisymmetric sausage waves that travel up and down the fibril. 

The paper is organized as follows: in Section [2] we present our scattering formalism. In 
Section |3J we calculate the scattering matrices. Section @] describes the absorption achieved 
and the properties of the mode mixing across mode order that occurs, while Section |5]presents 
the resulting phase shifts. Finally, in Section [6] we discuss our findings and conclusions, 
including what our results for a single tube adumbrate for plage. 



2. Scattering Formalism 

We consider a compact magnetic fibril that is embedded within a field-free atmosphere. 
The fibril is axisymmetric, straight and aligned with gravity. Under such conditions the 
acoustic wavefield in the surrounding, nonmagnetized atmosphere can be expressed as the 
sum of three components, 
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$(x,t) = $ inc (a;,t) + ^ 8ca (x,t) + $ jac (x,t) . (2.1) 

Here $ can represent any physical variable we might choose to consider; however, we will 
restrict our attention to a plane-parallel atmospheric model that is neutrally stable to 
convection — a good approximation within the sun's convection zone. For such an atmo- 
sphere, acoustic waves can be described using a displacement potential and we will let $ 
be this potential. The first term, $i nc , is the unperturbed p-mode wavefield that would 
exist in the absence of the magnetic flux concentration (i.e., the incident wavefield). The 
second term, $ sca represents a discrete set of far-field scattered waves in the form of propa- 
gating p modes. Finally, $j ac is the contr ibution due to the con tinuous spectrum of laterally 



evanescent, near- field, jacket modes (see lBogdan fc Calfylll995l ). 



We can express an arbitrary wavefield of incident p modes as a Fourier-Bessel expan- 
sion in cylindrical polar coordinates, x = (r,<j),z), where the origin is centered on the flux 
concentration, 



$ inc (u;;r,0,,2) = ^ ^ A ™(u) ^ J m (k n r) Q n (u; z) . (2.2) 

m=— oo n=0 

In this expression — and in the subsequent two — we have chosen to apply temporal Fourier 
transforms (t — > u>) and to work in frequency space. Our Fourier convention is such that 
our solutions have a temporal dependence of e~ luJt . The coefficients A mn (uj) are arbitrary 
complex amplitudes that characterize the incident wavefield, J m {k n r) is the Bessel function of 
the first kind, Q n (u; z) is the vertical eigenfunction for the nth-order p mode, and k n = k n (u) 
is the wavenumber eigenvalue for the nth-order p mode with frequency u>. 

Since we have assumed that the magnetic fibril is axisymmetric, incident waves of az- 
imuthal order m only scatter into waves with the same order. Therefore, for such an ax- 
isymmetric scatterer, the scattering only occurs across mode order n, and the explicit form 
of the scattered waves can be represented as follows: 

oo oo oo 

$ sca (u;;r,0,z) = ^ A mn (u) e im * £ S^ n '(u) Hg\k n ,r) Q n ,(u; z) , (2.3) 

m=— oo n=0 n'=0 

oo OO 

$ jac (u;;r,0,z) = ^ ^i„He' ra M dA T"~* a (uj) K m (Ar ) q(u, A; z) . (2.4) 

m=— oo n=0 ^ 

In these expressions Hm\k n /r) is the Hankel function of the first kind — corresponding to 
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an outward propagating cylindrical wave, and K m (Ar) is the modified Bessel function of 
the second kind, with A -1 being the lateral decay length for the jacket modes. There are 
two separate scattering matrices: S^ n ' is the far-field scattering matrix that represents 
scattering into outgoing, propagating p modes, whereas is the near-field matrix and 

represents the excitation of the acoustic jacket. In this formalism, mode mixing appears as 
nonzero values for the off-diagonal elements of the far-field scattering matrix (i.e., S^ n ' ^ 
for n 7^ n'). The functions q(u,A; z) describe the vertical behavior of each jacket mode and 
are oscillatory with depth. We note that the p mode eigenfunctions Q n (z) and the jacket 
modes g(A; z) form a complete, orthogonal set. We shall make use of this orthogonality 
property in the following section to compute the scattering matrices. Appendix |X] describes 
in detail the p mode eigenfunctions and jacket waves appropriate for the neutrally-stable, 
polytropic atmosphere. 

In general, the two scattering matrices, S^ n ' and T^ A , tell us everything we need to 
know about the wave-tube interaction, including the absorption coefficient a mn and phase 
shift A mn . For a particular mode (i.e., a particular m and n) the absorption coefficient is 
defined as one minus the ratio of the outgoing power to the ingoing power contained by that 
mode, 



a. 



10J 



A mn ^ (U) , 



Atl(uj) 



A 



(out) 



A m n (uj ) 



(2.5) 



The phase shift is defined in the usual way, as the phase difference between the ingoing and 
outgoing components, 



^ A-mn \UJ) J 

If we consider only a single incident wave, we can relate the absorption coefficient 
and phase shift to the scattering matrices in a straight forward manner. Noting that an 
incident wave component is comprised of outgoing and ingoing parts of equal amplitude, 
J m = (^Hm + HmJ /2, the outgoing and ingoing amplitudes, A^ out ^ and A^ m \ are then 
given by, 



= ^4^1 (2-7) 
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A^\u) = ^^[l + 25r»]. (2.8) 

If we now insert these expressions for the ingoing and outgoing power into equations (12. 5p 
and (12. 6p . we find 



a, 



10J 



1 + 2ST n (u) 



-4Re{5r n }"4 



Ql. 



(2.9) 
(2.10) 



A m » = arg{l + 2Sr»} , (2.11) 

2Im{gr n } ' 
,l + 2Re{^"}. 



tan" 1 ( 2Im {^r"} (2 ^ 



For realistic wave fields (e.g., as sampled by helioseismic observations) these simple 
formulae are invalid. Firstly, multiple incident waves of all mode orders n exist simultane- 
ously. Secondly, mode mixing (the scattering from one mode order n to a different mode 
order n') will result in the outgoing component of any mode order being comprised of the 
scattered waves from all incident mode orders. One goal of this paper is to assess the im- 
portance of mode mixing in the energy budget of an incident wave. Fortunately this can be 
accomplished by examining each incident mode in isolation, and computing the off-diagonal 
elements of the scattering matrix. We will leave for a subsequent study the effects that a 
complex, multi-component incident wave field has on measured absorption coefficients and 
phase shifts. Further, we will only consider mode mixing across mode orders. We pre- 
sume that the scatterer is axisymmetric; therefore, scattering from one azimuthal order m 
to another m! does not occur. Finally, for simplicity we have chosen to examine only the 
axisymmetric component (m = 0) of the incident wavefield, leaving higher-order waves to 
subsequent studies. The index m will be dropped from most expressions from here forward. 



3. Calculating the Scattering Matrices 



Scattering matrices have been co mputed previously fo r axisymmet ric magnetized tubes 



within unstratified atmospheres (e.g.. iKeppens et al.lll994l ; iTirryl |2000| ) . However, calculat- 



ing the scattering matrices for a magnetic tube in a gravitationally stratified atmosphere is 
formidable and can only be accomplished in general through numerical wave-field simula- 
tions. However, if we assume that the magnetic flux tube is thin and untwisted, substantial 



-7- 



progress can be made without resorting to computer simulation. By thin, we mean that 
all characteristic scale lengths — such as the wavelength of the external acoustic oscillations 
and the density scale height — are much larger than the radius of the tube R(z). Such a 
thin flux tube, although in hydrostatic balance, is unable to support internal lateral forces; 
hence, the total pressure will be uniform across the tube and continuous with the external 
value. This mandates that the total pressure has the same scale height inside and outside 
the tube. If we further assume that the temperature is the same within the tube as outside, 
the plasma parameter 0, defined as the ratio of the gas pressure to the magnetic pressure, is 
constant with height inside the tube. For photospheric and deeper layers, this assumption 
that the tube is in thermal equilibrium with its surroundings is consistent with observations 
( ILagg et al.ll2010l ). Thin tubes with constant /3 are theoretically attractive, as they permit 
analytic treatment of the waves that propagate along them and prove to be weak scatterers 
for which asymptotic analysis may be used. 

The scattering matrices are determined by the application of two physical matching 
conditions. At the interface between the magnetic tube and the surrounding field-free at- 
mosphere, both the total pressure perturbation and the normal component of the fluid dis- 
placement must be continuous. In general, thin flux tubes permit two types of magnetosonic 
waves, sausage and kink waves. The sausage waves are axisymmetric (m = 0), longitudinal, 
slow surface waves with pressure as the primary restoring force. The kink waves are sinuous, 
transverse, fast waves with m = ±1 and a restoring force composed of magnetic tension and 
buoyancy. Since we will only be considering axisymmetric incident waves, the solution within 
the magnetic tube will be composed of only sausage waves. We mention for completeness 
that magnetic tubes also permit torsional Alfven waves which are also axisymmetric oscilla- 
tions. However, acoustic-gravity waves within a neutrally-stable atmosphere are irrotational 
and are unable to couple to torsional waves. 

The propagation of sausage wave s on magnetic fibrils is reasonably well understood (e.g., 
Defouwlll976t iRoberts fc Webblll9781 ). Such waves have the property that the total pressure 
is nearly constant across the cross section of the tube and matches the overpressure in the 
field-free surrounding media. Furthermore, the normal component of the fluid displacement 
at the boundary of the tube is a small quantity of the same size as the radius of the tube 
(which is presumed small by the thin tube approximation). The smallness of the normal 
displacement enables the solution inside the tube to be obtained to leading order with only 
the enforcement of press ure continuity; the cont i nuity of the displacement appears only in 
higher-order corrections ( jSpruit fc Zweibellll979t ISpruitlll982[ lAndries fc Callyil201lh . Here 
we explicitly consider these higher order terms in order to compute the scattered wave field. 



Outside the flux tube, a small argument expansion of the Bessel functions appearing 
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in equation (I2.2p reveals that, to lowest order in the radius of the flux tube R, the pres- 
sure perturbation of the incident wave is a zero-order quantity 0(R°), whereas the normal 
displacement is first order O(R). Thus, the ordering of terms for the incident waves is the 
same as the sausage waves. However, similar expansions of equations (I2.3P and (12 .4p re- 
veal that the far-field and near-field scattered waves have a pressure perturbation whose 
leading-order terms are proportional to |S"|lni? or |T|lni2, respectively, where \S\ and \T\ 
indicate the leading-order behavior of the scattering matrices. The normal displacements 
for these waves are proportional to [S'l-R -1 and |T|i? _1 . The size of each of these terms for 
each wave component is summarized in Table 1. Careful scrutiny reveals that the manner in 
which a balance can be successfully achieved, where both continuity conditions are satisfied, 
is if the scattering matrices are to leading order proportional to R 2 . With this behavior, 
pressure continuity is maintained by a balance between the sausage wave and the incident 
wave; the scattered waves play no role to leading order. Continuity of the normal displace- 
ment involves a balance between the sausage wave and all three external wave components 
(incident, far-field scattered and near- field scattered). 

Table 1: LEADING ORDER OF EACH WAVE COMPONENT 



Wave Component 


Pressure Perturbation 


Normal Displacement 


External Incident Wave 


R° 


R 


External Scattered Wave 


\S\bxR 


\S\R-* 


External Jacket Wave 


\T\lnR 


\T\R- 1 


Internal Sausage Wave 


R° 


R 



By these arguments, the sausage waves are excited by the incident wave through pres- 
sure continuity and the scattered waves (and the scattering matrices) are generated by the 
sausage waves through continuity of the normal component of the displacement. The normal 
displacement outside the tube N ext is a linear combination of wave components, 

°° A. (to) 

N ext (u; z) 

,71 

(u; z) + iVj (u;z)} , (3.1) 

n=0 Z ° 

where iVl nC) „ is the normal displacement for the incident wave of order n and zq is a scale 
depth (chosen to be the depth of the photosphere in our field-free atmospheric model — see 
Appendix |A|). The displacements N sca ^ n and A^ aCin are the far- field and near- field scattered 
waves generated by incident mode order n. All displacements are evaluated at the flux- 
tube interface r = R(z) and any given displacement can be computed from the associated 
displacement potential through the relation, 
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N = n ■ V$ = — - — — 



I H if 



T -1 



(3.2) 



where n is the unit vector normal to the surface. The leading order behavior of each of these 
terms is found by small argument expansion of the Bessel functions in equations fl2.2l) - fl2.4p . 



N- mCin (uj; z) 

N sca , )n (uj] z) 
N iaCt n(u; z) 



-R{z)4 



zl 2i 



-2 Qn ^ z) + R^^z- {uJ;z) 



^ ' n'=0 



R(z) 



(3.3) 
(3.4) 
(3.5) 



Each order of incident wave excites a separate sausage wave on the magnetic fibril, 



N l] (u;z)=J2^t 1 N ln (u;;z). (3.6) 

n=0 

Appendix [B] provides details on how to compute each N\\ tn (u; z). Continuity of the normal 
displacement requires that iV|| >n = A^ inCi „ + A scain + iVj aCjri . If we insert equations H3.3I) — 
(13.51) into this relationship and multiply by R/zq, we obtain an equation for the "mismatch" 
between the incident and tube waves, 



d n (u;z) = [N- mCjn (uj; z) - Ni n (u; z)] , (3.7) 



Y J Sr n \u)Q n ,{u-z)+ / dA Tq~^ a (u) q(u), A; z) . (3.8) 



Figure [T] illustrates the two displacement components for an incident P2 mode, the compo- 
nents for the resulting sausage wave, and the mismatch between them d n . 

We can now solve for the two scattering matrices separately by invoking the orthogonal- 
ity of the eigenfunctions and the jacket waves. In Appendix |A] we discuss that the weighting 
function for the orthogonality integrals is proportional to the mass density. If we make use 
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of a dimensionless depth s = —z/zq, the density within a neutrally stable atmosphere is 
a power law p = poS a , where p is the photospheric density and the polytropic index a is 
related to the ratio of specific heats 7 through a = (7 — l) -1 . Multiply equation (13. 7p by 
s a and the complex conjugate of either a p-mode eigenfunction Q or a jacket wave q and 
integrate over all depths to obtain, 

S^ n '(u) = —J dss a Q n ,(u ] s)d n (u ] s), (3.9) 

/oo 
ds s a q*(u,A;s) d n (u;s) . (3.10) 

The limits of integration correspond to the photosphere (s = 1) and infinitely deep into the 
solar interior. 

Each element of the scattering matrices can therefore be thought of as the projection 
of the mismatch between the incident and scattered waves onto each p-mode and jacket- 
mode eigenfunction. The matrix associated with the far-field scattering, Sq~ Hi (u) can be 
easily obtained by direct numerical integration of equation ( 13. 91) . We have used a Bulirsch- 
Stoer numerical integrator to perform these integrations and the results are discussed in the 
following sections. The integral that must be computed to obtain Tq^ a (oj) is fundamentally 
more difficult to numerically compute since jacket wavef unctions are not confined to an 
acoustic cavity and remain propagating to any depth. Therefore, while the magnitude of the 
integrand does vanish as s — >■ 00, it does so slowly as a power law. This is unlike the integrand 
of Equation (13.91) which behaves like the product of a power law times a rapidly decaying 
exponential as s — > 00. Therefore, we delay the calculation of the near- field scattering 
matrix to a subsequent time. Fortunately, the far-field and near-field scattering matrices are 
independent of each other due to the orthogonality of the p modes and the jacket modes; 
thus, using this technique to calculate the far-field scattering matrix avoids the necessity of 
computing the acoustic jacket simultaneously. 

Under many reasonable boundary conditions, one can demonstrate that the far-field 
scattering matrix Sq ~ >n ' is symmetric under exchange of n and n' (see Appendix [D] for a 
derivation). Here we apply a stress-free condition at the photosphere (for both the p modes 
and the sausage waves) and finite energy (for the p modes) and radiation (for the sausage 
waves) boundary conditions deep in the atmosphere (as s — > 00). These boundary conditions 
are among those for which the scattering matrix is symmetric. This means that we need 
not calculate all of the matrix elements, only a little more than half need to be computed 
directly with the remaining obtained by symmetry. 



4. Absorption and Mode Mixing 



The absorption coefficient is often considered to be a measurement of energy loss or 
dissipation. When mode mixing occurs, this simple scenario needs to be modified as in 
addition to true absorption (or removal of energy from the acoustic wavefield), energy can 
be shuffled between modes by mode mixing. In this section we calculate the energy budget 
for a single incident wave and account for all the ways that the incident wave energy can be 
reprocessed by the scatterer. When the scatterer is an axisymmetric thin flux tube, there 
are only three ways in which energy can be extracted. (1) The single incident wave can 
excite tube waves which subsequently propagate out of the acoustic cavity. This mechanism 
is a true absorption process whereby acoustic energy is removed from the system. (2) The 
incident wave can generate an acoustic jacket that surrounds the flux tube. This mechanism 
is only significant for initial value problems and not for a spectral treatment as performed 
here. (3) The incident wave is scattered into outgoing waves of a different mode order. This 
last effect is mode mixing which is not true absorption as the acoustic energy remains in the 
cavity, albeit at a different wavenumber. We will assess the relative importance of the first 
and third mechanisms and leave the second to subsequent work. 



4.1. Energy Budget of a Single Incident Wave 

Consider a single incident wave of order n, with a unit amplitude (A mn = 1). This 
incident wave generates a set of outgoing waves with mode orders n'. These outgoing waves 
can be written as follows, where we have included both the far-field scattered waves and the 
outgoing component of the incident wave, 



n'=0 

The term in square brackets is the amplitude of the outgoing wave of mode order n', and 
is composed of a contribution from the incident wave (5 nn >/2) and from a scattered wave 
(Sq^™'). If the vertical eigenfunctions Q n are normalized according to equation f ]A6j) . each 
outgoing wave carries a lateral energy flux F n i with the following form, 



H£>(k n ,r) Q n .{u;z) 



(4.1) 



1 s 

Fn< = 2 Po u z o 



°nn' T aOq 



(4.2) 



Note, the energy flux only depends on the mode order n' (and the azimuthal order m) 
through the amplitude. Therefore, the fraction e nn i of the ingoing wave's energy that is 
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re-emitted in mode order n' has similar dependence on the scattering matrix as exhibited by 
the absorption coefficient, 



+ 2S^ n 



(4.3) 



For this example with only a single incident mode order, the absorption coefficient is by 
definition 1 minus the diagonal elements of this fractional energy flux, 



a n = 1 - e r 



|l + 2^ n | 2 



(4.4) 



-4Re{5^ n } -4|5 n ^ n | 



Of course this equation is identical to equation ( 12 .101) . By energy conservation, we can 
decompose this absorption coefficient into two parts, a true absorption from the excitation 
and escape of sausage waves W n and a redistribution term arising from mode mixing M n , 



a n = W n + M n . 



(4.5) 



The contribution from mode mixing M n can be directly written down using the off-diagonal 
elements of the fractional energy loss £ nn /, 



(4.6) 



Therefore, by combining equations (14.41) and (14.61) . the loss of energy due to the excitation 
of tube waves must have the following expression, 



-4Re{S^ n } -4^ \s't 



n— >n 




(4.7) 



n'=0 



Before we can use any of these expressions, we need to consider the order of each of their 
terms. Since we have computed the scattering matrix to only leading order in the radius 
of the tube, O (-R 2 ), for self-consistency we must drop the last term in equations (I4.4p and 
(HID, 



a n (u) = -4Re{^ n } + 0(R 4 ) , 
W n (cu) = -4Re {S^ n } + 0(i? 4 ) . 



(4.8) 
(4.9) 
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To lowest order, the absorption coefficient and the excitation rate of tube waves are identical. 
Therefore, our calculation of the absorption coefficient is only compatible with tube-wave 
generation. The terms due to mode mixing appear at higher orders which are not represented 
self-consistently. Fortunately, we can evaluate the mode mixing terms separately using the 
far- field scattering matrix through equations (14. 3[) and (14. 6p . 

Figure [2] illustrates the absorption coefficient (solid curves) evaluated through the use 
of equation (I4.8p . The fraction of the incident energy, M n , mixed into all other modes is 
represented with the dotted curves. In Figure [3J the energy e nn i that is scattered into each 
individual mode is shown. From this figure, one can see that the amount of energy lost to 
mode mixing decreases as the order of either the scattered mode n' or the incident mode 
n increases. The symmetry of the scattering matrix leads to the symmetry of e nn i, which 
is apparent in Figure |3fr as the equal height of pairs of symbols of different color and in 
Figure [3fr as the overlying of curves. 

For low-frequencies, the mode-mixing losses are rather inconsequential compared to the 
sausage-wave losses. However, for higher frequencies the amount of energy lost to mode 
mixing can be a relative high fraction of the total energy loss, perhaps reaching as high as a 
third for some parameters. Furthermore, for the / mode at high frequencies the excitation 
of sausage waves becomes anomolously inefficient and the mode-mixing terms can be the 
dominant source of energy loss. A more in depth discussion of tube wave excitation by the 
/ mode appears later in §6.11 



The computation of the scattering matrix also allows us to estimate phase shifts A n 
through equation (I2.12p . These phase shifts are illustrated in Figure HI For p modes, both 
positive and negative phase shifts are possible, although all shifts become positive for suffi- 
ciently high frequency — which, depending on the value of f3 may be beyond the physically 
relevant band of frequencies. The / mode is fundamentally different. Not only is the phase 
shift for the / mode negative for all frequencies, but the phase shift is large, reaching values 
as high as 15° for high frequencies and for tubes with low values of the plasma parameter /3. 

The phases of all the scattered waves are correlated in a precise manner with the phase 
of the incident wave. The difference in phase A0 nn / between each individual scattered wave 
and the incident wave can also be derived from the scattering matrix, 



5. 



Phase Shifts and Phase Correlations 
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These relative wave phases are uniformly small (see Figure [5]). The difference is positive for 
n + n' even and negative for n + n' odd. The symmetry of the scattering matrix to exchange 
of n and n' appears in Figure as the overlying of curves and the equal height of symbols. 
In general, scattering to and from the / mode results in larger phase differences, exceeding 
scattering between p modes by an order of magnitude. 

6. Discussion 

We have calculated by semi-analytic means the scattering matrix for axisymmetric 
acoustic waves encountering a vertically-aligned, magnetic tube. We have assumed that 
the tube is slender and invoked the thin flux-tube approximation to compute the scattering 
matrix to leading order in the tube's radius R for mode orders up to n = 10. In the following 
two sections we discuss the implications of the results of our calculation, both in terms of 
the energy budget of the incident and scattered waves and in terms of the phases of the 
outgoing waves. 

6.1. Absorption and Energy 

When encountering a thin flux tube of the type discussed here, a p mode of a particular 
radial order n may lose energy in either of two ways. The incident acoustic wave can excite 
MHD tube waves on the magnetic fibril that freely propagate along the tube carrying energy 
out of the acoustic cavity. This represents a true absorption as acoustic energy is removed 
from the wavefield. The tube may also scatter the incident p mode into outgoing waves of 
different radial order n' . Mode mixing of this sort does not represent a true absorption as 
the energy is merely redistributed between radial orders. Despite this fact, mode mixing can 
modify the observed absorption coefficient, as the definition of this coefficient does not take 
power redistribution into account. These two mechanisms are both inherently weak for thin 
tubes. The thin tube approximation mandates that the excitation of tube waves generates an 
absorption coefficient that scales linearly with the magnitude of the scattering matrix, which 
in turn scales like the square of the tube's radius R 2 , a small quantity by definition. The 
energy losses due to mode mixing are even smaller, scaling as the square of the magnitude of 
the scattering matrix, or as R . Therefore, mode mixing can only be important energetically 
when other sources of absorption (i.e., tube wave excitation) are peculiarly small. Over the 
physically relevant range of frequencies (< 6 mHz), our calculations reveal that mode-mixing 
dominates the absorption coefficient only for high-frequency / modes interacting with low-/? 
tubes. 
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Any model with a reflective upper boundary can have special frequencies for which the 
excitation of tube waves is unusually tiny. If the upper boundary reflects all wave energy, the 
excitation can formally be zero at these frequencies. Such "excitation nulls" arise because of 
total destructive interference. A localized driver, such as an incident p mode, excites waves 
that propagate both up and down the tube with equal energy. At these special frequencies, 
the upward propagating wave reflects of f the upper surface with just the phas e needed to 



cancel the downward propagating wave (IJain et al.ll2009l ; ICrouch fc Callyill999l ). In such a 



circumstance the downward energy flux vanishes and no energy is carried away from the 
acoustic cavity by tube waves. Near such an excitation null, the energy losses arising from 
mode mixing could be the dominant source of energy loss — not because mode mixing is large, 
but because the other source of absorption is small. 

For the stress-free upper boundary condition applied here and for reasonable atmo- 
spheric and tube parameters, excitation nulls do not appear within the physically relevant 
frequency band (u/2tt < 6 mHz). However, formally for the / mode such a null exists in 
the limit of infinite frequency and it has extremely broad influence, causing the absorption 
coefficient to asymptote to zero as the frequency becomes large. The frequency regime for 
which this asymptotic behavior is valid depends on /3, with lower /3 pushing the behavior to 
lower frequencies. Physically, such behavior exists because the / mode is a surface wave that 
becomes increasingly confined to the photospheric layers as the frequency increases. The / 
mode's evanescence length k~ x can be derived directly from its dispersion relation uj 2 = gk, 
revealing that the vertical extent of the / mode scales as u~ 2 . The wavelength L of the 
sausage wave at the photospheric surface can be derived from Equation (IB12I) . 



(7 - 1)3 " 



L (1 + 7/9/2K 



1/2 

(6.1) 



As the frequency increases, L shrinks, scaling as uj^ 1 ; however, it decreases more slowly 
than the scale length of the / mode decreases (which scales as u~ 2 ). Therefore, as the 
frequency increases, compared to the tube wave the driver looks more and more like a delta 
function located at the surface. A delta-function driver generates waves that propagate in 
each direction with identical amplitude and phase. If one of those components subsequently 
reflects off a surface and passes back through the driving layer, there is the potential for 
destructive interference (as discussed previously). Destructive interference occurs if the phase 
of the reflected wave advances by an odd integer factor of n. This phase advance can come 
from a phase change on reflection or from propagation of the wave to the reflecting layer and 
back. For the / mode, one can show that the phase change on reflection asymptotes to —tt 
as the frequency becomes large (see Appendix [U]). Therefore, the / mode, in the limit of 
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infinite frequency, satisfies the necessary resonance condition because the distance between 
the driving layer and the photospheric reflection layer vanishes. 

The p modes do not evince such high-frequency behavior for the simple reason that 
at high frequency, the depth of the cavity becomes roughly constant with frequency. The 
wavenumber eigenvalue k n (u) scales linearly with the frequency, thereby causing the hori- 
zontal phase speed u/k and the lower turning point of the p mode to approach constants. 
Since, the depth of the cavity and hence the vertical extent of the driving region does not 
change as the frequency changes, the wavelength of the sausage waves becomes much shorter 
than that of the p modes. Thus, for incident p modes, the driving region remains vertically 
distributed and the arguments that we made previously for the / mode are not valid. 

Figure M illustrates this behavior for (3 = 0.1. While the figure extends to frequencies 
as large as 30 mHz, one must bear in mind that frequencies greater than 6 mHz are greater 
than the sun's photospheric acoustic cutoff; thus, such waves are no longer trapped and 
the construction of our solar model is suspect. Further, one must be careful to not over- 
interpret the results at such extreme frequencies as the perturbation scheme has become 
invalid (see Figure E]). Despite this, for the / mode, trapped waves of low frequency (where 
the perturbation analysis remains valid) can still sense the presence of the excitation null at 
infinity, causing a high-frequency fall-off of the absorption. The frequency of waves affected 
depends on the plasma parameter (3 of the tube because the wavelength of the tube wave 
depends on (3. We can estimate the lowest frequencies that should be affected by requiring 
that the wavelength of the sausage waves L be two or three times larger than / mode's 
evanescence length k^ 1 . If we require that the ratio of these length scales satisies kL > 2, we 
discern that for (3 = 0.1 we would expect the high frequency regime to exist for frequencies 
above 3 mHz, whereas for /3 = 10 the regime begins around 10 mHz. 

This discussion only holds for absorption and mode mixing resulting from the excitation 
of the slow sausage wave. Similar calculations performed with fast kink oscillations reveal 
that excitation nulls can appear for all incident mode o r ders a t frequencies less than 6 mHz 



(IJain et al.l 120091 ; ICrouch fc Callyl Il999l ) ; also see iTirryl (120001 ) for a similar phenomenon in 



an unstratified model. Further, the calcu lations of the scatter ing; matrix elements arising 
from coupling to kink oscillations (S^ n ) ( Hanasoge et al.ll2008l ) show significant oscillatory 



modulation in the magnitudes as a function of frequency, whereas our scattering calculation 
for the axisymmetric sausage wave (iSg ~ m ) show smooth frequency variation. In fact one 
could fit the magnitude of the scattering matrix calculated here with a single power law in 
frequency that would be roughly valid over all but the highest frequencies. 

In a real plage, we do not expect such excitation nulls to have an easily identifiable effect 
other than a diminution of absorption at high frequencies for the / mode. Not only does 
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mode mixing tend to fill in the absorption profile, but the plage is likely to be comprised 
of a large number of tubes with a distribution of properties, including photospheric radius 
Rq and plasma parameter /3. Therefore, each tube in the collection will have suppressed 
absorption over a different ran ge of frequencies . The added affect of all of the tubes will be 



to mitigate the effect of nulls (IJain et al.l l201lh . In general, we expect most flux tubes in a 



plage to have a near equipartition of magnetic and gas pressure, thus tubes with (3 « 1.0 
should predominate. A quick examination of Figure |2] reveals that only the most extreme 
frequencies are likely to be affected. 

To summarize, we find that for thin tubes, mode mixing is a small effect compared to 
tube wave excitation. The scattering matrix itself has elements that are proportional to the 
cross-sectional area of the tube. Tube wave excitation extracts energy at a rate proportional 
to the scattering matrix (or proportional to the cross-sectional area of the tube) while mode 
mixing redistributes energy at a rate proportional to the square of the scattering matrix. 



6.2. Wave Phases 

Phase shifts are generally defined as the difference in phase between the outgoing wave 
in the presence of the scatterer and the phase that would have existed in the absence of the 
scatterer. A magnetic tube can generate a phase shift by three distinct mechanisms: (1) The 
ingoing acoustic wave reflects off the lateral surface of the tube, shortening the total path 
length traveled by the wave. (In the absence of the tube, the wave would travel all the way 
to the coordinate axis and back.) Since the wave travels a shorter distance, the presence of 
the tube causes a phase lag resulting in a negative phase shift or travel-time difference. (2) 
The wave is at least partially transmitted through the lateral surface and propagates at a 
different speed within the magnetic tube. If the phase speed increases, a negative phase shift 
is achieved. (3) Upon reflection off the lateral surface of the tube, the wave can undergo a 
phase change. Off course if a flux-tube boundary could be treated as a rigid, impenetrable 
surface, the phase change would be 180 degrees. But, in general a flux tube is elastic and 
the phase change may differ from 180 degrees and may be either positive or negative. 

For our thin flux-tube model, the only wave permitted within t he tube is a sausage wav e 



corresponding to a slow surface wave (see iRoberts fc Webblll978l ; lEdwin fc Roberts! Il983l ). 



For such a wave, the thin flux-tube approximation precludes lateral variation by implicitly 
assuming that the wave crossing time is much shorter than the wave period. Thus, for all 
tubes the crossing time for the tube should be ignored, resulting in a common negative phase 
shift, independent of /3. Obviously the phase shift whould be equal to the phase lag that 
would have occured in the unmagnetized model for a wave propagating a distance equal to 
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the tube diameter. This shows that for thin flux tubes, mechanisms 1 and 2 discussed in 
the previous paragraph are essentially the same. This is not true for thicker tubes for which 
the travel time across the tube is significant. A rough estimate of the expected phase shift 
would be —2kR , where i?o is the photospheric radius of the tube and we have assumed that 
the phase lag should be evaluated at the photospheric surface. This is clearly appropriate 
for the / mode as it is a surface wave confined to this layer. For p modes one should 
probably consider a depth weighted average of the tube radius. However, for our purposes, 
the previous expression will suffice. Estimates made this way for the / mode give a phase shift 
of Ao = —2u 2 Ro/g, which for waves at 3 mHz give a phase shift of -15°. This surprisingly 
large shift is about a factor of four larger than the phase shifts we obtain with our model (see 
HJ), indicating that a phase change upon reflection is likely to be at work as well. We point out 
that despite the large size of the estimated shift, it is quite similar in m agnitude and of the 



same sign as actually measured phase shifts for small magnetic features. iDuvall et all (120061 ) 
measured the average scattering from roughly 2500 small magnetic elements, and generated 
a typical travel time kernel from their measurements. One can obtain a typical phase shift 
of -14° by the follow i ng pr ocess. Multiply the largest value of the measured kernel in Figure 



6 from IDuvall et al.1 (120061 ) (roughly -20 s kG _1 Mm -1 ) by the average magnetic flux of the 
studied features (0.6 kG Mm 2 ) to obtain a travel time lag of -12 s. Convert this time lag 
into a phase shift by dividing by a wave period (300 s) and multiplying by 360 degrees. 

The phase change that occurs at reflection from the lateral boundary of the flux tube 
is difficult to estimate, as the reflection depends on both the incident wave and the driven 
sausage wave over a range of heights. We would like to receive guidance from calcula- 
tions performed for uns tratified, thin tubes, for which analytic solutions are possible (e.g., 



Edwin fc Roberts! 119831 ). However, the regime in frequency-wavenumber space for which 
slow surface waves exist and for which propagating external acoustic waves exist are dis- 
joint. Thus, for unstratified tubes a propagating acoustic wave cannot excite slow surface 
waves. The stratification in our model is what makes the coupling possible. We can safely 
say, however, that the phase shift caused by mechanisms 1 and 2 always result in negative 
shifts for thin tubes. Further, we would expect the phase change upon reflection (mechanism 
3) to increase for tubes of higher j3 since such tubes are more elastic. For the / mode this is 
indeed what we achieve; the phase shift is large and negative at all frequencies (probably as 
a result of a reduced path length) and becomes less negative as (3 increases (probably due 
to an increasing positive phase change at reflection). However, for p modes we obtain shifts 
of both signs; negative phase shifts at low frequency and positive shifts at high frequency. 
Clearly for the p modes a positive phase change upon reflection at the flux-tube boundary 
must occur for at least high frequencies. 
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A. The Field-Free Atmosphere and its Eigenfunctions 

The field-free atmosphere that surrounds the magnetic flux tube is treated as a plane- 
parallel, neutrally stable polytrope. This atmosphere has been utilized extensively in the 
past as its a good representation of a stellar convection zone and simplifies the wavefield 
by imposing a condition of irrotationality. Gravity is uniform, acting in the downwards 
direction g = —gz, with the height coordinate z increasing upward. The atmosphere is 
polytropic below the height z = —z , which corresponds to the model's photosphere. The 
gas pressure, mass density and sound speed vary with height as power laws. If we employ a 
nondimensional depth s = —z/zq, the background thermodynamic variables become, 
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The quantities p and P are the photospheric values of the mass density and gas pres- 
sure. The value of the polytropic index a is set such that the stratification is neutrally stable 
to convection; this requires a = 1/(7 — 1), where 7 is the ratio of specific heats. Above the 
photosphere z = —z , we assume the existence of a hot vacuum (p ext — > with T ext — > 00), 
with the property that the gas pressure (P ext ~ Pext^ext) is finite and continuous across the 
z = — z layer. 



Following iBogdan et al.l ( 119961 ) and iHindman fc Jainl (120081 ) , we specify the depth of 
the photosphere z and the photospheric density p ( an d therefore the referenc e pressure P n ) 
by ma tching our model photosphere to the r 50 oo = 1 level of a solar model by iMaltby et al. 
(119861 ). At this layer in the solar model g = 2.775 x 10 4 cm s -2 , p = 2.78 x 10~ 7 g cm -3 , 
and Pq = 1.21 x 10 5 g cm -1 s~ 2 . We adopt a polytropic index of a = 1.5, which is consistent 
with an adiabatic index of 7 = 5/3. 
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Since the buoyancy frequency in such an atmosphere is zero by definition, acoustic 
waves propagating within the atmosphere are irrotational. This permits such waves to be 
expressed using a displacement potential If we assume oscillatory time dependence e~ lU)t 
and harmonic functions with wavenumber k in the horizontal direction, the vertical portion 
of the acoustic eigenfunction satisfies the following Sturm-Liouville equation, 
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where A is a dimensionless horizontal wavenumber, v is a dimensionless frequency and Q 
represents either a discrete eigenfunction Q n {v; s) for real A = A n = 2k n z or a spectrum of 
jacket modes q(u, A; s) for continuous, imaginary values of A = — 2iAz . 

The general solutions to this equation are proportion al to Whittaker M and W functions 
( iWhittaker fc Watsonlll952l : lAbramowitz Stegunlll964h . The discrete / and p modes reside 
within an acoustic cavity and therefore become evanescent deep in the atmosphere. Thus, 
the W functions are relevant as they decay exponentially as s — > oo, 



Q n (u; s) = C n {v) s -(a^V2) W Rn ^{\ n s) 
Here C n {y) is a normalization constant and 



(A4) 
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The quantization of the horizontal wavenumber k n (or equivalently A n ) arises from the re- 
quirement that the Lagrangian pressure perturbation vanish at the model photosphere. 

The laterally evanescent jacket waves are surface waves (with a lateral decay length of 
A" 1 ) that propagate up and down the outside of the flux tube. They are not trapped within 
the acoustic cavity and can freely propagate to any depth. If we select only those jacket 
waves that satisfy the condition that the Lagrangian pressure perturbation vanishes at the 
model photosphere, the jacket waves are standing waves that can be represent ed as a sum 
of Whittaker M functions with imaginary argument (see lBogdan &: CaUylll995l ). 
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q(u, A; s) = C{v, A) s^ 1 /*) [r^(A) M Kfl (\s) + M K ^(\s)} , (A5) 
where C(u, A) is a normalization constant and 
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Unlike the / and p modes, for which the dimensionless horizontal wavenumber A is purely 
real and only takes on discrete values A n , for the jacket modes A is purely imaginary and 
continuous. Therefore, the argument and K-index of the Whittaker functions in equation ( 1A5I) 
are purely imaginary as well. 

Since both the discrete modes and the jacket modes satisfy Hermitian boundary con- 
ditions (vanishing Lagrangian pressure perturbation at both the photosphere s = 1 and at 
infinity s — > 00), we automatically know that the discrete, acoustic modes and the jacket 
modes are all mutually orthogonal with a weight function s a (which is proportional to the 
mass density), 



/oo 

ds s a Q n (u; s) Q n '{v] s) = 5 nn t , (A6) 

00 

ds s a Q n {v;s) q(v,A;s) = 0, (A7) 

ds s a q*(v,A;s) q(v,A';s) = 8 (A - A') . (A8) 



In writing these last equations, we have chosen C n (u) and C(v>,\) such that the discrete 
modes and jacket modes are orthonormal. 



B. Sausage Waves 

We have assumed the magnetic fibril is untwisted, straight, axisymmetric, vertically 
aligned and thin. Thin flux tubes are unable to support internal lateral structure and hence 
the total pressure is uniform across the tube and continuous with its external value. If we 
further assume that the tube is the same temperature as its surroundings, the total pressure 
must have the same scale height inside and outside the tube and the plasma parameter (3, 
defined as the ratio of the gas pressure to the magnetic pressure, is constant with height 
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inside the tube. Since we only consider tubes below the photosphere (i.e., z < —z ), we need 
not worry about the rapid flaring of the tubes into a magnetic canopy that occurs within 
the chromosphere. 

We ignore lateral variation of the magnetic field strength, and describe the tube's in- 
ternal gas pressure P(z), mass density p(z), sound speed c(z), and field strength B(z) by 
their axial values. These four quantities as well as the tube's radius R(z) can be described 
uniquely by the total magnetic flux contained by the tube and by the plasma (3, 
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Observations indicate that typical field strengths f or small phqtospheric flux tubes are 
betw een 1 and 2 kG within both quiet sun and plage (e.g.. lLagg et al.ll201ulMartmez Pillet et al 
19971 ). Such field strengths indicate rough equipartition between magnetic and gas pressure. 
For a flux tube with (3 = 1 embedded in the polytropic atmosphere described in Appendix [Aj 
the magnetic field strength at the photosphere is B = 1.2 kG. By inserting the external 
pressure P ext into the expression for the cross-sectional area of the tube flB5p . one obtains 
the explicit dependence of the tube radius on depth, 



R{s) = RqS 
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(B6) 



1 /2 

where Ro = (Q/nB ) 1 is the photospheric radius of the tube. For all illustrations we will 
assum e -Rn = 100 km, a fairly typical value for the size of magnetic elements within solar 
plage (ILagg et all 120101 ) . 



Thin flux tubes support a longitudinal magnetoacoustic wave known as the sausage 
wave. Here, these oscillations are driv en by the pressure per turbation imposed by exter- 



nal /- and p-mode oscillations (e.g., iRoberts fc Webbl 1 1 9 T8h . Using the formulation of 
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Hindman fc Jainl (120081 ) . to lowest order in the radius of the flux tube, the vertical displace- 



ment (\\ within the tube can be described by the following equation, 
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where ct is the cusp or tube speed which depends on both the sound speed c and Alfven 
speed Va within the tube, 



c 2 Vl T , a B 2 



c 2 + Vl' Va 4vrp 



The right-hand side of Equation (1B7[) is due to the forcing by the incident wave and is 
evaluated at the flux tube's axis, r = 0. This forcing arises from the enforcement of continuity 
of total pressure across the boundary of the flux tube. We are allowed to consider only 
the incident waves in this pressure balance because the far-field and near-field scattering 
amplitudes are smaller than the incident wave by a factor of R 2 In R. 

The forcing provided by each individual incident p mode can be obtained by evaluating 
the expression for the incident wavefield, Equation ( 12. 2p . at r = and inserting the result 
into the sausage wave equation. After Fourier transforming in time t, one obtains 

2 d 2 igc 2 T d 2 \ Pext 4 2 f , / dQ n (u-z) 

~ Y^d~z + U S C " = — Vf £ ^^tz— ■ (B8) 

Note that only the p modes with m = contribute to the forcing of the sausage waves. All 
other azimuthal components vanish in the limit r — > 0. 

We adopt a non-dimensional form for this equation by the following substitutions, 

z 2 + 7/3 7 (/3 + l)z/ 2 dQ n (z/;s) 

s = > e = o — ' -k( s ) = o 5 
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The solution of this equation is expressed using a Green's function, 
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In this equation, if)\\(s) and V^ 5 ) are the two solutions to the homogeneous equation, J n (s) 
and X n are interaction integrals over the driver, and Q n is a parameter that specifies the 
complex amplitude of the wave that reflects from t he upper surface and th us determines the 
boundary condition at the model photosphere (see lHindman &: Jainll2008l . for details), 
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If we apply a stress-free boundary condition on the sausage wave at the photosphere, 
the boundary condition parameter has the following form, 



. 7(1 + 1) 2 Q n n* 

tt n = z v — — — X n , (B15) 

7r rt rt 

H = uV~eH^ 1 (2uV~e) + (P + l)(fi + l)H^(2uV~e) , (B16) 

On = Qn{v,s)\ s=1 . (B17) 

Figure [TJi, b illustrates the resulting sausage wave as well as the driving p-mode eigenfunction 
when the incident wave is a P2 wave. 

We can derive an expression for the normal displacement inside the tube using the 
definition for the divergence of the displacement vector £ in cylindrical coordinates, 

„ dh\\ h\\ dC\\ ,„ , 

V • = + ^ + , (B18) 
or r oz 

where h\\ is the radial displacement. In this equation we have utilized the axisymmetry of 
the sausage mode. This can be combined with the Taylor series expansion of the radial 
displacement, 
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which we evaluate at r = -R(-z) to obtain 
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If we now use Equation A9 from iBogdan et al.l (119961 ) , which is a statement about pressure 
equilibration of the tube with its external environment, 
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and combine equations (13. 2p and (1B20j) . we obtain after some manipulation 
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This displacement can be decomposed into seperate components driven by each incident 
mode, by utilizing Equations (12. 2p and (IBlOj) . 
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C. High-Frequency Asymptotics for Incident / Modes 

The eigenfunction for the / mode is simple enough that we can perform high-frequency 
asymptotics to derive approximate analytic solutions for the interaction integral Iq and all 
quantities derived from it. The / mode eigenfrequency is always u 2 = gk independent of the 
details of the stratification. The associated eigenfunction is just an exponential that grows 
with height Qo(u; z) = exp(kz). In our dimensionless variables, these relations take the form 
v 2 = a\/2 and Qo(v; s) = exp(— v 2 sjd). The / mode's evanescence length is given by ajv 2 
and as the frequency becomes large, the evanescence length becomes tiny, scaling as ~ v~ 2 . 
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The / mode's eigenfunction is, therefore, progressively more confined to the photospheric 
surface as the frequency becomes large. The tube wave on the other hand has a wavelength 
that diminishes much more slowly. From equation (IB12I) the wavelength of the sausage wave 
can be seen to scale as ~ v~ l . Thus, in the interaction integral X — see equation (IB14[) . the 
integrand is only significant near the lower limit of integration (near the photosphere) and 
the sausage wave can be considered long wavelength. 



C.l. The Interaction Integral 

After inserting the /-mode eigenfunction and the homogeneous sausage wave solution- 
equation (1B12|) — into equation (1B14|) . we obtain 
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Since the exponential term in the integrand has significant magnitude only very close to 
the lower limit of integration and the sausage wave is long wavelength compared to the 
exponential scale length, we may expand the function F(s) in a Taylor series about s = 1, 

F(s) = F(l) + F'(l)(s - 1) + *-£l{8 - l) 2 + • • • (C3) 

If we define 
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and perform the necessary derivatives, the following expression is obtained, 



F(s) = + - ue 1/2 h^ +1 ) (s - 1) 



+ - 



2 2 

-/i — jj, — v e 



h,. 



(/i - l)ue 1/2 h 



1) 2 + -- 



(C5) 



After inserting this expansion into the integral appearing in equation (IClft . each term 
in the expansion can be integrated separately by analytic means, producing the following 
expansion for the interaction integral, 
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x = 7(/3 + 1) ,V- 2 /- 



hp + ^ (/i^ - ^ 1/2 /^ M +i) - + 0(i/ 7/2 ) 



(C6) 



The order of the missing terms was deduced by noting that for large frequency the large- 
argument expansion of the Hankel function in equation (1C4I) is appropriate. Therefore, 
hp ~ v~ x l 2 and the leading order frequency dependence for Zq is 0(v 3 / 2 e~ u2 / a ). To lowest 
order, the interaction integral should decay rapidly as the frequency increases. 



C.2. Amplitude of the Downward Propagating Sausage Wave 

The complex amplitude of the wave that reflects from the photosphere is given by 
the boundary condition parameter Q. For the stress-free upper boundary with an /-mode 
incident wave, the boundary condition parameter Qq takes on the following form, 



.7(0 + 1) 



-u 2 /a 



If we now combine equations (1C6|) and (1C7j) we find asymptotically, 



(C7) 



7(0 + i; 



n 



2i 

7i 



+a(fi + l)((3 + l)u- 1 e 1 / 2 h^hp +1 
+a 2 p- l e^% +1 hp + O (1s- 3 ) 



-1 l/2r* l 

a\w e 1 h^ +1 hp 



(C? 



From equation (IB 1 1 j) we can easily see that well below the driving layer, as s — > 00, 



the amplitude of the downward propagating sausage wave becomes proportional to fi +Xq 



Oddly, if we add expression (IC8I) to the complex conjugate of equation (1Q6|) and collect 
terms of like order, all of the terms cancel to the order calculated. Therefore, for an incident 
/ mode, in the limit of infinite frequency, the reflected wave and the directly generated wave 
undergo complete destructive interference and infinite frequency should be considered an 
"excitation null" . 

In this limit, the driver is essentially a delta function located at the surface. Therefore, 
the phase difference between the reflected wave and the directly generated downward wave 
can only accrue from a phase change upon reflection off of the photosphere. For destructive 
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interference, this phase change must be ±7r. The actual phase change turns out to be 
AO = —it, and this can be verified asymptotically. The phase change can be expressed using 
the ratio of the complex amplitudes of the upward propagating wave and the reflected wave. 
These can be obtained by evaluating equation (IB11I) at s = 1, 

A " = arg {§}' (C9) 

Since the frequency is large, the Hankel functions represented by can be expressed using 
their large argument expansions, 

AO « arg ( ^L e -^ 1/2 e <n+W>> \ . ( C 10) 



If we now compute the ratio Iq/Qq using these same large argument expansions we find, 

^0 _ e 4ii/e 1 /2 e _i( M+ 3/2) 7 r _ {Cll) 

Combine these last two expressions to find that AO = —it to leading order. 



D. Symmetry of the Far-Field Scattering Matrix 

In this appendix we will demonstrate that under several commonly used boundary 
conditions the far-field scattering matrix, Sfi -m , is symmetric to interchange of n and n! 
(S$-> n ' = S^ n ). We start with equation 1ET5]) . 

/oo 
dss a Q n ,(s)d n (s) (Dl) 

where the frequency dependence has been suppressed. First we need an expression for the 
mismatch, d n , which is supplied by equation (13. 7J) , 

dn(s) = ^ [N iuc>n (s) - N ln (s)] . (D2) 



Using equations (Q, flEfi} and flB24j) this expression can be reduced to a sum of two 



terms: a term that depends solely on the incident p-mode eigenfunction and a term that 
depends on the excited sausage waves, 
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d n {s) = ^s~ a { s» 



s ds 4 



2s 



7/3 d 
2sz ds 



(*"C||,») 



(D3) 



The terms in the square brackets can be modified using the Sturm-Liouville differential 
equation that the p-mode eigenfunctions satisfy, equation ( 1A1I) . 



dJs) 



Rl s -a 



2zl 



(D4) 



Insert this expression into equation (IDlj) to obtain, 



on—¥n 
D 



ff {°n'n + ^ r ds Qn>{s)i [^M] 



2ez J t ~ u " y ds 
where the values D n > n are the elements of a differential operator, D 



(D5) 



DQ(s) 

Dn'n 



^ WW] - ^J^Qi*) , 



0?S (5n'( S ) D< 5"(s) • 



(D6) 
(D7) 



Since we have applied Hermitian boundary conditions to the p modes, one can easily 
verify that the operator D is also Hermitian with respect to the same set of functions; thus, 
D n / n = D nn '. So, we now only need consider the second term in the braces in equation ( 1D5I) : 
Is this term symmetric upon exchange of n and n'l 

The longitudinal displacement C||,n can easily be expressed using a Green's function, 
G{s,s'), 



C||,r 



^0 



/ ds'G(s,s')f n (s') 



(D8) 



where 



G(s,s') 



(D9) 
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/ n \Ms)^!l(s') H S'<S, 

I ip*As) ip\\(s ') if s' > s, 

/.(.) = _20L+JL1^. (mi) 

2 s as 

This formulation is different from the one used in Appendix B because we will exploit 
symmetries possessed by the Green's function and the present formulation makes those 
symmetries explicit. In particular, note that the function \l/(s, s') is symmetric upon the 
interchange of s and s' . 

In expression ( 1D8|) . Q n is a complex constant that determines the boundary condition 
applied to the sausage wave at the surface of the truncated polytrope. For a stress-free 
surface the appropriate value of fl n is given by, 

^tW + i)^.^ (D12) 

71 ft ft 



and, alternatively, for an upper atmosphere with minimal reflectivity ( IHindman fc Jainll2008l ) 
the correct value of the boundary condition parameter Q n is 

n n = - l{(5 + l) v 2 KQ n . (D13) 

If one inserts equation flD8j) into expression flD5j) . integrates by parts, and uses the fact 
that the p-mode eigenfunctions and the tube wave solutions vanish at infinity (s — > oo), one 
can obtain the following expression after some manipulation, 

Oq —— — —U n i n + — j [Wn'n + ±S n 'n) , 

4 Zq L'OEZq Zq 

where 



HV „ s -AiiAs r r*>« , (Di5) 



i ji 



Bn ' n = 7 (^+ X i) t / 2 " Sn ' + • (D16) 

It is trivial to demonstrate that W n > n is symmetric. Since the variables s and s' in the 
integrals are dummy variables, we can exchange them (s — > s' and s' — > s) and then exchange 
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the order of integration. Further, the function ^(s, s') is symmetric to interchange of its two 
arguments. Therefore, after performing these exchanges we find the forestated symmetry, 
W n , n = W nn .. 

Thus, the symmetry of the far-field scattering matrix hinges on the properties of the 
matrix B n / n , which contains boundary terms arising from the integration by parts as well 
as a contribution that comes from the homogeneous solution appearing in the first term 
in equation (1D8|) . Under some common boundary conditions for the sausage waves, this 
matrix is symmetric, but in general it is not. One can readily verify by directly inserting 
expressions flD12j) or (ID13j) into equation (ID16j) . and using the identity for the Wronskian of 
two Hankel functions, 



W{HS?{x),Hg\x)} = hS +1 {x)H^{x)-H^\x)H^ +x {x) = -— , (D17) 

that both of the boundary conditions mentioned previously (stress free and minimum reflec- 
tivity) produce a symmetric matrix, B n / n = B nn i. Therefore, at least for these two boundary 
conditions, the far-field scattering matrix is also symmetric. 

The crucial properties that were required for this symmetry are (1) the p modes satisfy 
Hermitian boundary conditions, (2) the sausage wave vanishes deep in the atmosphere (s — > 
oo), and (3) the sausage wave satisfies a boundary condition at the surface of the truncated 
polytrope (s = 1) that causes the boundary term matrix B n i n to be symmetric (note, not 
necessarily vanishing). This proof is very reminiscent of the calculation that is performed 
to demonstrate that a differential operator is Hermitian; however, in this case the operator 
in question is actually an integro-differential operator (because of the integration over the 
Green's function). 
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Fig. 1. — (a) The vertical displacement £i nc (black curve) and the normal displacement to the 
flux tube's surface -/Vj nc (red curve) for an incident p2 mode with a frequency of 3 mHz. (6) The 
vertical and normal displacements (black and red curves, respectively) for the sausage wave that 
is driven by the incident p2 mode. The magnetic flux tube has a plasma parameter of (3 = 1. (c) 
The resulting mismatch in the normal displacement between the p mode and the sausage wave. In 
all panels the solid and dashed curves respectively represent the real and imaginary parts (when 
present). 
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Fig. 2. — Fractional energy lost from the incident wave as a function of frequency for 100 km radius 
flux tubes with the three indicated values of the plasma parameter (3. The colors indicate the order 
of the incident waves (/: black, p\: red, p2- green, p^: blue, p^. aqua, p§: violet). The solid curves 
are the absorption coefficient a n computed with the scattering matrix using equation (|4.8|) . The 
dotted curves correspond to the energy lost from the incident mode order due to mode mixing into 
all other mode orders, M n , equation (|4.6p . The energy loss arising from mode mixing is generally 
small compared to the tube wave excitation, although at high frequencies (> 5 mHz) the mode 
mixing can climb to one-third of the total energy loss. Further, for the / mode in the limit of 
large frequency, the absorption asymptotes to zero, leading to the property that the mode mixing 
can dominate the energy loss. The transition into this regime occurs when the wavelength of the 
sausage wave becomes much larger than the evanescence length of the / mode; thus, the transition 
frequency increases as j3 increases. 
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Scattered Mode Order - n ' Frequency - co/2ji {mHz} 

Fig. 3. — Fractional energy lost from the incident mode due to mode mixing induced by a magnetic 
tube with plasma parameter /3 = 1 and a radius of 100 km. (a) Energy lost to each mode order 
of scattered waves for an incident wave of 3 mHz in frequency. The labels and the different colors 
correspond to the order of the incident wave, using the same color scheme as in Figure 2. (b) 
Energy lost to mode mixing as a function of frequency. The colors indicate the incident mode 
order (/ through ^2)1 while the linestyles correspond to different orders of the scattered modes (/: 
solid, p\\ dots, p2- short dashes, p^: dot-dashed, p^: three dots and a dash, p^. long dashes). In 
both panels careful examination verifies that the modulus of the scattering matrix is symmetric, 
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Fig. 4. — The phase shift as a function of frequency for flux tubes with a radius of 100 km and 
with the three indicated values of f3. The different colors represent different incident mode orders 
(/: black, p\\ red, p2' green, p^: blue, p^. aqua, p$: violet). The bottom panels show the same 
quantity, but with a differently scaled ordinate in order to better illustrate the / mode. Both 
positive and negative phase shifts are possible for p modes and are universally small. Whereas for 
the / mode, the shifts are negative and substantial, reaching values as large as -15°. 
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Scattered Mode Order - n ' Frequency - co/2ji {mHz} 

Fig. 5. — The phase of the scattered waves measured relative to the phase of the incident wave, 
(a) Absolute value of the phase of each mode order of scattered waves for an incident wave of 3 
mHz in frequency. The x symbols correspond to positive phases and □ symbols to negative phases. 
The different colors (and labels) indicate the order of the incident wave, (b) Relative phase of each 
scattered wave as a function of frequency. The colors indicate the incident mode order (/ through 
P2), while the linestyles correspond to different orders of the scattered modes as in Figure 3 (/: 
solid, pi: dots, p2'- short dashes, p^: dot-dashed, ^3: three dots and a dash, p^: long dashes). 
In both panels careful examination verifies that the phase of the scattering matrix is symmetric, 
arg{5^"'} = arg{S ™'->"}- 
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Fig. 6. — The high-frequency behavior of the absorption profiles. The thick solid curves are the 
absorption coefficient calculated via the far- field scattering matrix with equation (|4,8p . with the 
different colors indicating the order of the incident mode. The / mode's high-frequency behavior 
is fundamentally different from the p modes because it is a surface wave. At high frequency, 
driving from the / mode becomes spatially localized at the surface, wheras the driving from p 
modes remains spatially distributed. The thin, dashed curves show the modulus of the scattering 
matrix |iS r g~ >n (a;)|. For the perturbation scheme used here, the scattering matrix must be a small 
quantity. Clearly, for the / mode, this weak-scattering approximation is invalid for frequencies 
exceeding roughly 6 mHz (the shaded region). We show such large frequencies merely to illustrate 
the mathematical behavior at high frequencies. 




